A general model of watershed extraction and representation using globally optimal flow paths and up-slope contributing areas

Many modern hydrological models require data inputs provided by automated digital terrain analysis functions incorporated into GIS. These inputs include Ž elds representing surface  ow directions, up-slope contributing areas, and sub-catchment partitions. Existing raster-based terrain analysis tools, including both those in oŒ-the-shelf GIS packages and those in the recent literature, were designed to work with digital elevation data in mountainous topography. For highly variable topography, which may include large  ood plains, lakes, wetlands, and other relatively  at areas, existing tools cannot accommodate the variable signal-to-noise in the source elevation data without signiŽ cant human intervention to handle special cases. A general model for calculating  ow directions, up-slope contributing areas, and sub-catchment partitions that automatically adapts to the variable information content of grid-based elevation data sets is presentedhere. The model uses a combinationof breadth-Ž rst search and global optimization to extract the maximum amount of signal from any location within the data. The model is demonstrated to work well in handling topography dominated by large  ood plains, lakes and other  at areas without the need for a large number of empirical rules. An important contribution of the approach is the handling of explicit hydrologic features, whichmakes the spatial representation closely related to hydrological processes. The results have important implications for developing hydrological models that are tractable in large, heterogeneous watersheds using moderate resolution data.